Quantum trajectory approach to the geometric phase: open bipartite systems 
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Through the quantum trajectory approach, we calculate the geometric phase acquired by a bipar- 
tite system subjected to decoherence. The subsystems that compose the bipartite system interact 
with each other, and then are entangled in the evolution. The geometric phase due to the quantum 
jump for both the bipartite system and its subsystems are calculated and analyzed. As an example, 
we present two coupled spin-^ particles to detail the calculations. 
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Consider a quantum system that depends on some ex- 
ternal parameters X . We are interested in the evolution 
of its quantum states when the parameters X change 
slowly along a closed path. For an eigenstate, such 
an adiabatic evolution accumulates a geometric phase, 
known as the Berry phase Q], which reflects the sys- 
tem geometry with the parameter space X. Although 
the geometric phase has been rigorously formulated for 
the general case of non-adiabatic, non-cyclic, and non- 
unitary evolution of pure states, the importance of 
geometric phase in realistic systems, for instance, in the 
context of geometric quantum computing ^ l^i 
has motivated recent interesting research into geomet- 
ric phases for mixed state s B and open systems 
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Quantum trajectory (quantum jump) analyses have 
been applied to certain physical systems, which show how 
the geometric phase for a closed system can be modified 
under open system dynamics [l4l llq . The trajectory 
analyses can be applied for any systems evolving under 
a Markovian master equation or, equivalently, under any 
trace preserving completely positive maps, and it proves 
to be particularly suitable for the geometric phase be- 
cause in each particular trajectory the quantum state of 
system remains pure. However, the analyses presented in 
Ref.'T^ are only for a whole/single system, and those in 
Ref. 16] are for conditional phase gate (it is not a geo- 
metric phase) as well as for a single-qubit geometric phase 
gate based on numerical simulations. Then the problem 
of geometric phase in coupled open bipartite systems re- 
mains untouched. 

It is more important, from the perspective of quan- 
tum computing, to study the geometric phase in bipartite 
systems, since almost all systems employed to perform a 
quantum gate are composite, i.e., it at least consists of 
two subsystems with direct couplings or coupled through 
a third party. This together with the above motivation 
stimulate the interest in study of the geometric phase in 
open composite systems. 

This paper contains the following interesting advances 
in this field. First, we present a completely general calcu- 
lation for open bipartite systems with only one subsystem 
subjected to decoherence, showing the effect of decoher- 
ence on the Berry phase of the bipartite system. Second, 



having calculated the geometric phase for the composed 
subsystems, we show that quantum jumps occurring in 
one subsystem make no contribution to the geometric 
phase for another. Third, we identify what kind of jump 
operators does not change the geometric phase for both 
the bipartite system and its subsystems. Although the 
situation of only one subsystem subjected to decoherence 
seems less realistic, it leads us to see how decoherence in 
one subsystem affects the geometric phase of the whole 
system and of another subsystem, moreover the repre- 
sentation for this simple situation can be easily extended 
to the case of both subsystems subjected to decoherence. 

We start with the most general autonomous differential 
equation for the state of an open system in the Lindblad 
form ,2^] 

ih^^p=[H{X),p]+C{p), (1) 

where H{X) — H^{X) is the composite system Hamilto- 
nian depending on external parameters X. £(/9) repre- 
sents the Liouvillian, which has the general form l2,l ] 
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{TiTkP + pTin-2TkpTi}. (2) 



[H{X),p] in Eq.Q generates the coherent part of the 
evolution, while C{p) represents the effect of reservoir 
on the dynamics of the system, the action of each F^ 
amounts to a different decoherence process. Suppose that 
we monitor the system and do not detect any decay, the 
geometric phase for the no-jump trajectory of the master 
equation in the continuous limit is given by E] 



-dt-arg{«(r)|V/(0))}, (3) 



where i-^\'tp^{t)) — H\'ip^{t)), and H stands for the non- 



Hermitian effective Hamiltonian 



H = H{x)-'-J2riT, 



(4) 
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In the adiabatic limit, the geometric phase for a cyclic 
evolution among path |(/)„(X)) yields 12^] 
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with |<I>„(X)) and |0„(X)) satisfying 

H\4>n)^K\<t>n), and #t|$^j^A:j$„), (6) 

the parameter /argument X is omitted here and hereafter 
where it could not make confusion. Now, we generahze 
the formulation presented in for single systems to 
bipartite systems. When the bipartite system undergoes 
an adiabatic evolution along path the reduced den- 
sity matrix for one subsystem, say a (similarly for 6), is 
given by 

p^i,^Tn\cp.n){'^n\. (7) 

Having written state in the form of Schmidt decom- 
position 

J 

the reduced density matrix p" can be expressed as 

where \E^)a comes from the Schmidt decomposition for 
\^n). i.e., 

I'f") = E ® (10) 
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For the simple case where each pair of and are 
X-indcpcndent, the Berry phase of the bipartite system 
reduces to 

llab = ME V/P^(/^)*(71J + lnb,j)}. (11) 
j 

With 7l^, ^ il{E,{X)\V^e,{X))JX, and ^ 
i S^{Fj{X)\V^fj{X))bdX. It is easy to see that 7° ,^^, 
generally is not a weighted sum over the one particle 
geometric phases ^na.j ^^"^ Inbp since y/p^ji^/Pf)* are 
complex for open systems in general. This is different 
from the Berry phase in closed bipartite system 24] . 

Now, we turn to study the effect of the quantum jump. 
Suppose that the decoherence is only caused by a lo- 
cal reservoir, this indicates that each TkpT\, generates a 
quantum jump within one of the subsystem in the trajec- 
tory. Without loss of generality we assume here that the 
jumps occur only within subsystem a, i.e., all com- 
mute with any operator from subsystem b. In the quan- 
tum trajectory analyses, the dynamics of the bipartite 
system is approximated by dividing the total evolution 
time T into a sequence of discrete intervals 5t = T/N. 
The time evolution of the density matrix in each interval 
takes the form 16] pab{t + St) = J2k=o ''^kPab{t)w\, where 
wq = 1 — iH5t and — T}^\/5i (fc = 1, ...,n). If there is 
only one jump characterized by T in the trajectory at an 
arbitrary time ti, which occurs in a time much shorter 



than any other characteristic time of the system. Then 
the reduced density matrix pa after the jump reads, 

pl(ii) = Tn[Tm,))m,)\T^] = Tpa{h)T\ (12) 

where Pa{ti) represents the reduced density matrix of 
subsystem a before the jump. Then the phase associated 
with the occurrence of a jump at time ti is given by p5| 

7r"^ = ^,g{Y^p,{t,){a,{h)\T\a,{h))} 

j 

= arg{Tr[p„(ti)r]}. (13) 

Here Pa{t) = Pj 1*^^ (^)) ("^j (01 was used in the expres- 
sion, this is the geometric phase of subsystem a acquired 
in the jump, which is obviously non-zero. But the sub- 
system b acquires zero geometric phase associated to the 
jump, this can be understood as follows. The total phase 
shift due to the jump is 

iir' = ^rgmh)\rm,))}, (w) 

where \tp{ti)) denotes the state of the bipartite system at 
the time of jump. Have writing \tp{t)) into the Schmidt 
decomposition, \tp{t)) = J^j \/Pj\Pj{t))a'^\Tj{t))b, we ob- 
tain 

j 

= arg{Tr(p,(Or)}, (15) 

it is exactly the phase shift acquired by subsystem a as- 
sociated to the jump. Thus for a bipartite system, the 
subsystem that has no change under the action of the 
jump operators acquires zero geometric phase associated 
with the quantum jump, even if the coupling between the 
subsystems arc not zero. This result sharply depends on 
the assumption that the jump does not need time, i.e., 
it happens immediately and lasts no time. The situation 
changes when we lift the limitation/assumption on the 
jump, the jump in one subsystem would transfer a geo- 
metric phase to another due to couplings between them. 

To be specific, we apply this general representation to 
two coupled spin-i particles, in which one of the spin-i 
particle is driven by rotating magnetic fields and sub- 
jected to decoherence. We calculate and analyze the ef- 
fect of decoherence on the Berry phase of the bipartite 
system as well as the geometric phase of the subsystems. 
Let us start with the Hamiltonian that describes two cou- 
pled spin-i particles in time-dependent magnetic fields, 

H = ^aaa-B{t) + J{a+<j+ + h.c.), (16) 

where aj = (crj ,a^,aj), are the pauli operators 
for subsystem j{j — a,b) and cr+ ~ (l/2)(crj -|- icr]'). 

We will choose B{t) — Bon{t) with the unit vector 
•h — (sin^ cos sin0 sin0, cos^) and have assumed that 
only the subsystem a is driven by the external field. The 
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classical field B{t) acts as an external control parameter, 
as its direction and magnitude can be experimentally al- 
tered. J stands for the constant of coupling between 
the two spin-i particles. This coupling is not a typi- 
cal spin-spin coupling, but rather a toy model describing 
a double spin flip; nevertheless, the presentation in this 
paper may be generalized to the system of nuclear mag- 
netic resonance (NMR) where we can use Carbon-13 
labelled chloroform in acetone as the sample, in which 
the single ^^C nucleus and the nucleus play the role 
of the two spin-i particles. The constant of spin-spin 
coupling Jcrf erf in this case is J ~ (27r)214.5Hz, and we 
may control the rescaled coupling constant g = 2 J/aBg 
by changing the magnitude of the external magnetic 
field. We would like to address that the interaction 
between the two spin-i particles in our model is not 
a typical spin-spin coupling as that in NMR. So, we 
have to make a mapping when we employ the presenta- 
tion in NMR system and when all subsystem are driven 
by classical fields. The Liouvillian which describes the 
decoherence effect in subsystem a may have the form 
^P) = -t{<^^o-a P+Pcr+cr--2cr-pcr+}. The correspond- 



ing non-Hermitian Hamiltonian reads H = H — ^na^a~ , 
dissipation would give rise to modify the eigenvalues and 
the corresponding eigenvectors that are given by 



and 
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aj{(l},d,9)\eg) + bj{(j),e,g)\e 



Cj(0,6',5)lff5> + dj{(l),e,g)\ge)], 



(17) 
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aj{(j),0,g) 
dj{(j),0,g) 



9, g) — Ej — cost 



J4> 



bji(t>,0,g) = - 



g(cos 9 — Ej — |k) sin 9 
sin^ 9 ~ {cose ~^K)^ + Ef ' 
cos9 + Ej — 



sin0 



(18) 



El = ^lsm'9+icos9--K.)^ + -^ + -^^g^+4sm'9 = -E2, 



Es = ^/sin^ 61 -I- (cos6l- ^k)2 -f ^ - 1^52 ^_ 4 gij^2^ ^ 



(19) 



The eigenvectors and the corresponding eigenvalues of 
have the same form as those of H in Ea. (|18ll9|l . 
but {—in) should be replaced by in. We will use nota- 
tions of Aj,Bj,Cj, and Dj {j — 1,2,3,4) corresponding 
to aj , bj , Cj , and dj in Eq. (|18|l as the coefficients that 
appear in the instantaneous eigenstates of H'' . In fact, 
Aj{iK) — aj{—iK), Bj{iK) — bj{—iK) and so on. The 
numerical results for the Berry phase with no-jump were 
presented in figure 1, where we plot 7° ^jf, (n = 1, 2, 3, 4) 
as a function of the spontaneous rate k and the azimuthal 
9. In contrast with the case of k = [23|, there are jumps 
at 9 = t: /2 with k > kq depending on the path the sys- 
tem follows. The jumps appearing in figure 1 may be 
understood as follows. The evolution of the state along 
the no-jump trajectory represented by one of equation 
(I17|l can be mapped on the Bloch sphere with sponta- 
neous decay rate ^ k. The evolution is then a smooth 
spiral converging to the lower state, thus the geomet- 
ric phase increases due to the spontaneous decay when 
the initial state falls onto the upper semi-sphere; while 
the phase decrease when it initially is in the lower semi- 
sphere. Therefore, the geometric phase has a jump at the 
crossover point 9 = 7r/2. This was schematically shown 
in figure 2. In the case where only one quantum jump 
(described by operator a~) occurs at any time ti, the 



geometric phase shift due to the jump is given by 

7Ca7 = arg{a„C: + fe„i?:}, {n = 1, 2, 3, 4), (20) 

selected numerical results for 7,^™^ were illustrated in 
figure 3. The quantum jump occurring at a time of </> = tt 
was assumed for this plot. Clearly, singularity appears at 
9 — 7r/2, where the subsystem a experiences a crossover 
from the upper half of the Bloch sphere to the lower one. 

In conclusion, we have investigated the geometric 
phase in open bipartite systems. This study is of rele- 
vance to the geometric quantum computing, where ge- 
ometric phase may be used to perform quantum infor- 
mation processing. To our best knowledge, this issue 
remains unaddressed by the trajectory approach, in par- 
ticular for coupled open systems. The results show that 
there is a singularity in the dependance of the geomet- 
ric phase on the azimuthal angle with a specific k > kq, 
where kq depends on the inter-subsystem coupling con- 
stant g. The phase shift due to the quantum jump was 
also calculated, it has shown that the phase shift depends 
strongly on the direction to which the spin-^ particle 
points. The jump occurs ai 9 = 7r/2 can be interpreted 
in picture of the Bloch sphere, in which the decaying of 
the subsystem a results in a smooth spiral converging 
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FIG. 1: The Berry phase 7„,ab for the bipartite system in 
a no-jump trajectory as a function of the decoherence rate k 
and the azimuthal angle 9. (l)-(4) in the figure correspond to 
difi'erent instantaneous eigenstates j = 1,2,3,4. 



to the ground state, so when the state falls in the up- 
per semi-sphere, the decay increases the geometric phase, 
but it lowers the geometric phase when the state on an- 
other semi-sphere, this leads to the jump in the phase at 
9 = 7r/2. And it is interesting to note that there is no 
jump when the decay rate k < kq, the critical value of 
the decay, this is due to that the system finishes a cyclic 
evolution before the system decay. 



FIG. 2: A schematic illustration for the effect of spontaneous 
decay. The solid angle increases due to the decay when the 
state is on the upper sphere, while it decreases for that on the 
lower sphere. 
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FIG. 3: Geometric phase acquired by the bipartite system (or 
subsystem a) in a single quantum jump among the adiabatic 
path I'&j). The quantum jump was assumed to occur at the 
time when (j) = it. 



